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1. Introduction 

The summation of logarithms of \jx in deep inelastic structure functions at small 
values of Bjorkcn x leads to the Balitskii-Fadin-Kuraev-Lipatov (BFKL) equation 
[1, 2], which in the leading approximation sums terms of order [ctg ln(l/a;)]". Recently 
the next-to-leading terms have also been computed [3, 4]. 

Since the dynamics of the small-x region is supposed to be different from that 
at higher x in several respects, it is important to make and test predictions of a 
wide range of observables in this region. In recent papers [5, 6] predictions were 
presented for the rates of emission of fixed numbers of 'resolved' final-state gluons, 
together with any number of unresolvable ones. Here 'resolved' means having a 
transverse momentum larger than some fixed value ^i^. The predictions were valid 
in the double-logarithmic (DL) approximation, i.e. retaining only terms of the form 
[tts ln(l/a;) ln((5^/Ai|)]". In this approximation, each resolved gluon can be equated 
to a single jet, since to resolve it into more than one jet would cost extra powers of 
CKg with no corresponding powers of ln(l/2;). 

The present paper extends the work of ref . [6] to include terms with fewer powers 
of ln((5^///|), i.e. those of the form [cKg ln(l/a;)]"' [ln((5^///|)]"* where < m < n. 
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which wc refer to as single-logarithmic (SL) corrections. The fact that we still demand 
a factor of ln(l/a;) with each power of as means that the identification of resolved 
gluons with jets^ remains valid. 

There are two alternative methods for the calculation of final-state properties at 
small x: the original multi-Regge BFKL method and the CCFM [7, 8, 9] approach, 
which takes account of the coherence of soft gluon emission. It has been shown, 
at the DL level in refs. [5, 6] and now at the SL level [10], that the two methods 
are equivalent for the observables considered here. We therefore adopt the BFKL 
approach, which is calculationally simpler. 

The paper is organized as follows. In sect. 2 we recall the BFKL formalism and 
the predicted behaviour of the gluon structure function at small x. In sect. 3 we 
first compute the single-jet rate to SL accuracy and show that the only modification 
to the DL result comes from the SL corrections to the anomalous dimension. Next, 
in subsect. 3.2, we calculate the SL corrections to the two-jet rate, first in the form 
of a numerical integral and then as a perturbation series. In subsect. 3.3 we apply 
the same methods to the three-jet rate. In sect. 4 wc derive the SL perturbative 
expansion of the generating function for the multi-jet rates, and use this to obtain 
the corresponding expansions for the mean jet multiplicity and its dispersion. Our 
conclusions arc presented in sect. 5. A useful class of integrals is evaluated to SL 
accuracy in the appendix. 



2. BFKL formalism 



We start from the unintegrated structure function of a single gluon, /(x, /c^, /x^), 
which in the exclusive form of the BFKL approach satisfies the equation 

fix, k\ = 8{l-x) 5\k) + a, j^^ J A(^, k\ /.^) f{x/z, \q + k\\ /x^) . (2.1) 

Here a.s — 3Q;s/7r, k is the (2- vector) transverse momentum of the gluon probed in 
the deep inelastic scattering, q is that of an emitted gluon, is the azimuthal angle 
giving the direction of g, // is a coUinear cutoff and A is the Regge form factor 

/ I k'^\ 

A(z,/c^A^2) =exp(-asln-ln— j . (2.2) 

To carry out the z integration it is convenient to use a Mellin representation, 

Uk\ii^)= dxx-f{x,k\ii^), (2.3) 
Jo 

^When counting jets in deep inelastic lepton scattering, we always omit final-state hadrons that 
originate from the quark- antiquark pair which couples the gluon to the virtual photon. 
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with inverse 

= — / du x'^-' Uk\ ^i') , (2.4) 
zm Jc 

where the contour C is parallel to the imaginary axis and to the right of all singu- 
larities of the integrand. This gives 

Uk',^i') = 6\k) + H^{k',^i') f ^^f^i\q + k\',^i') (2.5) 

where 

To solve eq. (2.5) one can make use of the relation, derived in the appendix, 

(2.7) 

which is valid for /i^ < A;^ < to logarithmic accuracy, i.e. neglecting terms sup- 
pressed by powers of fi^ /k"^ or kP' /Q"^. Then for k"^ > ij? the solution is of the form 



where 



1 = H^ik', ^l')(ln-^-^ + 2J: C(2m + 1)(7 - 1)^"^ ) (2.9) 



and hence 

u = -as [27^ + ^^(7) + ^(1 - 7)] , (2.10) 

ip being the digamma function and 7e = —'0(1) the Euler constant. The solution of 
eq. (2.10) is 7 = 71^(0:3/0;), the Lipatov anomalous dimension: 



7.(as/a;) = ^ + 2C(3) p + 2C(5) p + 12[C(3)]^ p + . . . . (2.11) 
The integrated gluon structure function at scale Q'^ is then given by 

i^a;(g^//') = l+7^/ dA;v.(^^/i') = p . (2.12) 

Since wc arc interested in final states, we shall need to decompose the struc- 
ture function into terms corresponding to different numbers of emitted gluons. The 
contribution to f {x , k'^ , /s^) from emission of n gluons is obtained by iteration of 
eq. (2.1): 

f-\x,k',^') = f[ f ^^^a,A{z,,kln')S{x-x^)S\k-k^), (2.13) 
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where 

i i 

Xi^^zi, ki^-Y,Qi- (2-14) 

1=1 1=1 

The contribution to the structure function at scale Q is then obtained by integrating 
over all //^ < qf < Q^: 

fW(x, g^ ^^') = n r «sA(^„ n')S{x - x^) , (2.15) 

or in terms of the Mellin transform 

Fi-W,^^') = n /f ■ (2.16) 

3. Jet rates 
3.1 Single-jet rate 

Consider first the effect of requiring one emitted ghion, say the j'th, to have g| > //| 
while all the others have qf < This defines the contribution of one resolved gluon 
plus n — 1 unresolved, F*^"'^-*''*^: 

Notice that for i < j the contribution is identical to the {j — l)-gluon contribution 
to the structure function evaluated at = On the other hand for i > j we 
have /c? ~ g| > //|. As shown in the appendix, when g| > we can write for any 
function /, to logarithmic accuracy, 

/<M^,„...r)^^f (3.2) 

Thus the qi integrations for i > j become trivial and 

^ j=i ■'t'i ^3 

where we define^ 
Summing over all j and n gives the total one-jet contribution, 

Jix^^ qj 

= eMld^sMS] G^^\T) (3.5) 

^Note that we define 5 and T differently from Refs. [5, 6] (twice as large) in order to simplify 
expressions for the single-logarithmic terms. 



SH^iqlli^) 



n— j-t-i 



(3.3) 



4 



where 

G«(T)=ln(l + ^T) . (3.6) 

Notice that the colhnear-divergent part (the S'-dependence) factorizes, and the frac- 
tion of events with one jet is given by the cutoff-independent function 

Thus in the case of the single-jet rate, the only subleading logarithms are those 
generated by the presence of the full Lipatov anomalous dimension (2.11) in eq. (3.7). 
To obtain the jet cross section as a function of x, we note that eq. (3.5) implies 

that 

Fi''''\Q',^^l) = FM)Gj^\T), (3.8) 

where we have used the factorization property to replace the cutoff-dependent gluon 
structure function F^{ii\, /x^) by the measured structure function of the target hadron 
at scale //|, Fj^(/i|). It follows that the single-jet contribution as a function of x is 
given by the convolution 

F(^j^*)(x,g^/.^) = F{x,ixl)®G^'\x,T) ^ ['-F{z,i,l)^G^'\x/z,T) (3.9) 

Jx Z 

where the inverse Mellin transformation (2.4) apphed to eq. (3.6) gives 

G^''\x,T) = — / du;x-^-^G^^\T) = \ , , ^ . (3.10) 
^ ' ' 2mJc ^ ' a;ln(l/x) ^ ' 

3.2 Two-jet rate 

Now suppose we resolve two gluons j, / {j < f) with transverse momenta g|, g|/ > 
In place of eq. (3.5) we have 

(3.11) 

where, defining the dijet transverse momentum qj — qj + qji, 

KM'j,^l) = HM'j,^')(^+ t f[ [f^Y^HUklA . (3.12) 

When q^j > we can safely set kf = qj for i > j' , to obtain 

K^{qlfxl)=H^{qlfxl). (3.13) 

However, this cannot be correct for qj < fi'^, because K^{qj,fiD would then be 
infinite at qj — /x|exp(— tu/cis)- Since K^^i^q^j^ii'^ must be independent of the cutoff 
we can evaluate it for qj < /i| by setting /i^ = qj, which gives 

(2 \ 7l(^s/^) 
% . (3.14) 
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Thus we can define the continuous function 

K^, 111) = l^l)0{q] - lA) + -FM, QjWU - q'j) ■ 

The two-jet rate is then given by 

p(2jet)('/-)2 ,,2 ,,2\ 



where 



dqj dql dcpj 



(3.15) 



(3.16) 



(3.17) 

The integrations in eq. (3.17) can be performed numerically, without encounter- 
ing any non-integrable divergences or discontinuities in the integrand. The resulting 
two-jet rate is shown by the solid curve in fig. 1 for a relatively small value of ag/cj 
(0.2), and for a larger value (0.4) in fig. 2. The dashed curves show the result of 
using the DL prediction for G^^ in eq. (3.16), i.e. 



(3.18) 



where [6] 



^(2,DL) (J.) ^ A f 1 + ) + In ( 1 + ) - ( 1 + 



(3.19) 



We see that for a^/uj = 0.2 the single- logarithmic correction is small, while for the 
larger value it is substantial. 



0.15 



ijL'j"'^ 0.1 




0.05 



6 8 10 

T = ln(QVM|) 

Figure 1: Twojet rate for as/u = 0.2. Dashed: double-log approximation. Solid: with 
single-log corrections. 
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Figure 2: Two-jet rate for as/to = 0.4. Dashed: double-log approximation. Solid: with 
single-log corrections. 



Next we consider the perturbative expansion of the two-jet rate. We can use 
eq. (3.15) and the results in the appendix to obtain 



dq], d(j)j. 
Ml ql 27r 



■Ku,{\Qj + Qj'\ 



R ^J' 
OO 



+ 2 Y: (2m)! C(2m + 1) [H^iq], /.^)]^-+^(3.20) 



m=l 



Substituting in eq. (3.17) we find 



G(f)(T) = Gi2'°L)(T) + G(f'SL)(T) 



(3.21) 



where G^^'^^^ is the double-logarithmic result (3.19) and the single-logarithmic cor- 
rection is 



1 + 



-2m-l 



. (3.22) 



Note that the series in eq. (3.22) is strongly divergent for any value of as/uj. This 
is due to the singularity of H^{q^j^^'^) at gj = /x^ exp(— cj/ag) < The change in 
K^{qj, when < /x| (see eq. (3.15)) removes the singularity from the integrand 
in eq. (3.17) but this does not affect the region of convergence of the perturbation 
series. The situation is analogous to the way in which the running of the QCD 
coupling as{(f') at low values of (f produces infrared renormalons [11]: the Landau 
singularity in the perturbative expression for as{q^) leads to a factorial divergence 
of perturbative expansions with respect to q;s((5^), where is fixed and large, even 
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if wc remove the Landau singularity by making a non-perturbative modification of 
as(g^) at low [12]. 

In the case of eq. (3.20), the correction arising from the more careful treatment 
of the region qj < corresponding to the second term in eq. (3.15), would be 



I 

J u 



uj q'^ 




7L(as/a;) 



1 - ^j.(as/u;) 




(3.23) 



This does indeed contain a factorially divergent series, but, owing to the overall 
factor of 1/q], it does not contribute any logarithms of Q^/fi"^ upon substitution in 
eq. (3.17). Therefore a more accurate treatment of the region qj < although 
necessary to evaluate the integrals, does not affect the two-jet rate to SL precision. 

If we interpret the series in eq. (3.22) as an asymptotic expansion, then the 
partial sum truncated after the smallest term represents an estimate of the total SL 
correction, with an uncertainty of the order of the smallest term. This estimate is 
shown by the points in figs. 3 and 4, with the uncertainty represented by the error 
bars. We see that estimates from the series are of the same order of magnitude as the 
numerical results, but the discrepancy may be several times the expected uncertainty. 
At the smaller value of ots/oJ, the SL correction is relatively small (c.f. fig. 1) and the 
discrepancy is not so important. For the larger value oi a^luJ, the estimate is better 
than expected, but the correction and the uncertainty are both large. 



0.012 



0.006 




4 6 8 10 

T = \n{Q^/„l) 

Figure 3: Solid: single-log correction to two-jet rate for a.s/oj = 0.2. Points: estimate 
from asymptotic expansion; 'error bars' indicate the smallest term in the expansion. 
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0.02 



T = ln(gVM|) 

Figure 4: Solid: single-log correction to two-jet rate for as/uj = 0.4. Points: estimate 
from asymptotic expansion; 'error bars' indicate the smallest term in the expansion. 

To deduce the two-jet rate as a function of x, we can proceed as in eq. (3.9), 
writing 



where 



Zm Jc 

For the DL contribution we find from eq. (3.19) that 

1 



■e?(''°L)[asTln(l/a;) 



xln(l/x) 
where 

^(2,DL)[^j ^ ^^^^^ + + + 1 + [Ei(-^) + Inz + 7e - 1 - -z] 
£^1(2;) being the exponential integral function 

E,{z) = / dt— , 



(3.24) 
(3.25) 

(3.26) 
(3.27) 
(3.28) 



interpreted as a principal- value integral when 2; < 0. 

The divergence of the series for the SL correction in a;-space, eq. (3.22), is cured 
when one makes the inverse Mellin transformation to x-space, because the factorial 
coefiicients are cancelled: 



"i Ic 



1 

2m Jc 



, .^_,fasY^+'_as [«sln(l/x)]^^ 
'^''^ 'J ' X (2m)! 



UJ 



(3.29) 
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Thus the SL correction to the two-jet rate can be expressed in closed form as a 
function of x: 

G^''^^\x,T) = V7T7T^^''''^["Bln(l/a;)] (3.30) 
xln[l/x) 

where 

g{2,SL) ^ - ^) - In r(l + z)- . (3.31) 

Notice that the expression (3.30) is singular at Q;sln(l/x) = 1. Prom the viewpoint 
of the Mellin transformation (2.3), it is this singularity that produces the divergence 
of the series in cj-space. Conversely, use of the more correct expression (3.17) in 
a;-space, with the kernel function Ki_j given in eq. (3.15), should suffice to remove the 
singularity in a;-space. 

3.3 Three-jet rate 

The method used above for the two-jet rate can be extended, albeit laboriously, to 
higher jet multiplicities. In the three-jet case we have 

i^if /^l) = ^^''^^Ig^^^^j^'^ = eM-l.{^sMT] G(f)(T) (3.32) 

where 

j/"r qj qy qj" 

(3.33) 

One could in principle evaluate this expression numerically using cq. (3.15) for ii'^. 
Here we derive the perturbative expansion analogous to eq. (3.22). Introducing 
t — ln(g^//i,|) etc. for brevity, the results in the appendix give 



(3.34) 



L^{t)^H^{t)[hiH^{t)-hiH^{T)\+t[H^{t)f + 2 {2m')\C{2m' + l)[H^{t)\ 

m'=l 



G^^\T)^j dtH^it) I di'L,(max{i,i'}) + 2 5^C(2m + l)--^ 
JO [Jo ^^-^ at 

where we write H^{q^,ii'^) as H^{t) and 

oo 

2m' +2 
(3.35) 

Hence we obtain 

G^f ) (T) = Gif (T) + G(f 'S^) (T) (3.36) 
where is the double-logarithmic result [6] 

Cm^T) ^ 1 + ^t) + In^ (l + ^t) 

+ In f 1 + ^t) (i + ^tY- (l + (l + ^r) " (3.37) 

V UJ J \ UJ J UJ \ 2uj J \ UJ J 
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and 



oo /n, \ 2m+l 

^(3,SL)(j.) ^ -2$:(2m)!C(2m+l)p^ 



1 



2m + 1 I 

oo 

+ 4 E 



m=l 



a; 



-2m-2 



In 1 + + V'(2m + 1) + 7^ + 2 



1 + — T 



-2m-l 



+ 



(2m + 2m' + l)! ^.^ , ^, 



2m+2m'+2 



X 



-2m-2m'-2' 



(3.38) 



The expansion in eq. (3.38) is again strongly divergent for all values of cis/'^- As 
discussed in subsect. 3.2, it can still be interpreted as an asymptotic expansion and 
used as a guide to the order of magnitude of the SL correction. Furthermore, the 
divergence is cured upon inverting the Mellin transformation to obtain the jet rate 
in x-space, as long as asln(l/x) < 1. To extend the prediction to smaller values of 
X one would need to evaluate eq. (3.33) numerically using the full kernel function in 
eq. (3.15). 



4. Generating function for multi-jet rates 

For a general jet multiplicity r, we can write 



(4.1) 



u=0 



where the jet-rate generating function is given by 

R^{u, T) = exp[-7,(as/^)T] G^{u, T) (4.2) 

and 

00 

G^{u,T)^Y.^^G^;^\T). (4.3) 

r=0 

The function G[[)(T) was given in eqs. (3.6), (3.21) and (3.36) for r = 1,2 and 3, 
respectively. 

We can obtain the perturbative expansion of the function Gi_j{u,T) as follows. 
We first define the unintegrated function goj{u, t, T) such that 

G^{u,T) = 1+ r dtg^{u,t,T) . (4.4) 
^0 

Then using the results in the appendix we find that goj{u,t,T) satisfies the integro- 
differential equation 

1 + / dt' («, max{t, t'}, T) + 2 ^ C(2m + 1) 

■'^ m=l ^ 

(4.5) 



g^{u,t,T) = uH^{t) 
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Writing 

gUu,t,T)^J2('n{u,t,T)[^) , (4.6) 

n=0 ^ 

this implies that 

rT 

Cn+l{u,t,T) = u6n,0 ~ - u)tCn{u,t,T) + U dt' Cn{u,t' ,T) 



+2.EC(2m+l)^. (4.7) 

m=l 



Starting from cq = 0, this gives Cn iteratively as a polynomial in u, t and T, which 
can be substituted in eq. (4.4) to obtain the perturbative expansion of G^{u,T) to 
any desired order.^ The relation (3.29) can then be used to transform the result 
directly to x-space, giving 



F('■j«^)(x,g^/.^) = ^F{x,t,l)^—G{u,x,T) 



(4.8) 



u=0 



where 



with 



G(u, X, T) = (5(1 - x) + dt g(u, x, t, T) (4.9) 



as ^ [Qisln(l/x)]" 
g{u,x,t,T) = —}_^Cn+i{u,t,T) j , (4.10) 



which we beheve to be a convergent series as long as ag ln(l/x) < 1. 

4.1 Anomalous dimension 

Notice that ior u— 1 we have 

oo 

i?.(l,^) = Ei?(:J^*)(g^/.|) = l (4.11) 

r=0 

and therefore 

(7,(1,7) =exp[7,(asHr]. (4.12) 

To show that eq. (4.5) docs indeed lead to the Lipatov result (2.10) for the anomalous 
dimension, we note that when u = 1 the solution of eq. (4.5) is 

^,(l,t,r) = 7e^(^-*) 

G,(1,T) = 1+ r dtg^{l,t,T) = e^^ (4.13) 
Jo 



'The results agree with those given to fourth order in ref. [9]. 
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where 



7 = H^it) 



l + 7t + 2^C(2m+l)7 



2m+l 



m=l 



a; 



l + 2^C(2m+l)7 



2m+l 



m=l 



-7 



7 



- - V(l + 7) - - 7) 



(4.14) 



Rearranging terms, we obtain 7 = ^^{a^/uj) given by eqs. (2.10) and (2.11). 
4.2 Jet multiplicity moments 

We can compute the moments of the jet multiplicity distribution by successively 
differentiating the generating function at m = 1: 

Qs 



{r{r- l)...(r-s + l)) 



du- 



:Rco{u,T) 



(4.15) 



u=l 



In this way we obtain the perturbative expansion of the mean number of jets 



+4C(3)T^ (^)'-8C(5)t(|) 



+ 



(4.16) 



and the mean square fluctuation in this number, 



CO z 



- 2C(3)r 

40 



OJ 



+ 12C(3)T^ (^) - (8C(5)T - yC(3)T3) + . . . _ (4 ^7) 

It appears true to all orders to SL precision that, as in the DL approximation [6], 
the mean number of jets is a quadratic function of T and the mean square fluctuation 
is a cubic function of T. Thus the distribution of jet multiplicity at small x and large 
T is narrow, in the sense that its r.m.s. width increases less rapidly than its mean as 
T increases. 



5. Conclusions 

The calculation of jet rates at small x poses many interesting challenges and sheds 
new light on the novel dynamics of this kinematic region. In the present paper 
we have concentrated on those perturbative contributions which have a factor of 
ln(l/x) for each power of and are further enhanced by one or more powers of 
T — ln((5^///|), /^R being the minimum resolved jet transverse momentum. 
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For sufficiently large values of fi^ and ^ /^|, the resummation of such terms 
would seems to be a well-defined problem in perturbation theory. The results in 
sect. 4 do indeed specify all terms of the form {as/uj)'^T"^ with m > 0, for any jet 
multiplicity, cu being the moment variable in the Mellin transform. However, as we 
have seen explicity for the two- and three-jet rates (and we believe to be true more 
generally) , the single-logarithmic terms (those with < m < n) cannot be resummed 
directly since they form strongly divergent series. In ci;-space, the divergence is 
associated with kinematic regions in which the vector sums of transverse momenta 
of combinations of jets arc less than /Xr. A more careful treatment of such regions 
renders the jet rates well-defined as integrals. Furthermore, one obtains convergent 
series, within a limited range of x, after performing the inverse Mellin transformation 
to x-space. In the case of the two-jet rate, we were able to sum the resulting series 
explicitly, to obtain a closed-form expression valid in the region as\n{l/x) < 1. 

A number of interesting questions arise from our results. Clearly one would like 
to extend the resummation of jet rates to higher multiplicities and smaller values 
of X. This will require an x-space treatment of the difficult kinematic regions men- 
tioned above. One would also like to prove the conjectures in subsect. 4.2 about 
jet multiplicity moments to all orders, and preferably to resum them. Ultimately, 
next-to-leading terms in ln(l/a;) should also be included. Such terms will arise from 
next-to-leading corrections to the BFKL kernel and from the resolution of emitted 
gluons into two jets. 
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A. Useful integrals 

To evaluate integrals of the form 



/ 



f 27r 



/(k+^r) 



(A.l) 



to logarithmic accuracy, we may assume that /(A;^) has an expansion in powers of 
InA;^, and so it suffices to consider /(/c^) = In^/c^ for positive integer values of p. 
This can in turn be done by considering f{k^) — k'^" for small u and extracting the 
coefficient oiv^/pl. Now 



q^ 27r' 



i<j> 



21/ 



E 

m=0 



m 



fc2 



dq^ 



/ 2\rn 



dl_ 

jk2 g2 



A;2 



(A.2) 
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If k"^ < /i^ the first integral on the right-hand side is absent and the lower limit on 
the second becomes /x^; if k'^ > the second integral is absent and the upper limit 
on the first becomes Q^. Thus, neglecting power-suppressed terms, we have 



= k 



2v 



In 



1 

V 



for ji"^ < k"^ < Q"^, 



Here we have used the remarkable result, valid for Re > —1, 



(A.3) 



E 

m=l 



m 



1 



1 

— + 
m m — u. 



2 J2 C(2m + 1) = _27^-^(l + t/)-V;(l-t/) . (A.4) 



m=l 



Thus, neglecting power-suppressed terms, we find 



2 27r 



)2 J^2 



/(g^) for /c^ < /x^. 



+2 5^C(2m + l) A;^ 



m=l 
12^ 



_d_ 

dk' 



2m 



f{k^) for //2 < < 



Q 



In (^^1 /(A;') for k' > Q' 



(A.5) 
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